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Questions about black holes in quantum gravity generally presuppose the presence of a horizon. 
Recently Carlip has shown that enforcing an initial data surface to be a horizon leads to the correct 
form for the Bekenstein-Hawking entropy of the black hole. Requiring a horizon also constitutes 
fixed background geometry, which generically leads to non-conservation of the matter stress tensor 
at the horizon. In this work, I show that the generated matter energy flux for a Schwarzschild black 
hole is in agreement with the first law of black hole thermodynamics, 8nGAQ = kAA. 



INTRODUCTION 

In the search for a theory of quantum gravity, black 
holes play an important role. Since black holes behave 
as thermodynamic objects, with an entropy given by 
S = A/AHG 00, they presumably are connected to the 
underlying microstates of quantum gravity. Indeed, al- 
most every modern approach to quantum gravity makes 
attempts, with varying degrees of success, to show that 
the number of microstates for a black hole matches that 
predicted by black hole thermodynamics (see Q for a re- 
view). With a top-down approach this is feasible - one 
has the fundamental theory and can define the particular 
solutions that represent a black hole in the macroscopic 
limit. These solutions have a statistical mechanics de- 
termined by the fundamental theory, and the resulting 
macroscopic thermodynamics can be compared with the 
laws of black hole thermodynamics. 

Unfortunately, we don't have a complete theory of 
quantum gravity. It is therefore useful to analyze the 
statistical mechanics of black holes from a bottom up 
approach as well, in the hopes that information might 
be derived about the physics the statistical mechanics 
of the black hole microstates should follow. In doing so 
we immediately run into a puzzling problem. Black hole 
thermodynamics is well defined in the semiclassical ap- 
proximation: we know what spacetime represents a black 
hole and can investigate the behavior of quantum fields 
and small gravitational fluctuations on top of the back- 
ground spacetime. In a full theory of quantum gravity 
this splitting is unacceptable. There is no background 
spacetime available and metric fluctuations due to the 
uncertainty principle will prevent the localization of a 
horizon. Now our problem becomes clear: if we cannot 
even identify the states that represent a black hole, how 
can we possibly analyze their statistical mechanics? 

Two approaches have been developed for dealing with 
this issue. The first is the "horizon as boundary" ap- 
proach 4] which treats the black hole horizon as a bound- 
ary H. Imposing boundary conditions on the path inte- 
gral for the metric inside and outside the horizon such 
that H is a horizon yields the appropriate number of 
boundary states. Another, more recent, alternative is to 
impose constraints directly on the theory such that H 



is a horizon. In a recent work, Carlip has shown that 
imposition of "horizon constraints" yields the appropri- 
ate number of microstates for black hole entropy 0|. In 
this approach H is a spacelike initial data surface con- 
strained to be a stretched horizon. For a review of both 
these approaches see also 0] and references therein. 

The notion of a stretched horizon is important for this 
work, so I now take a little time to explain what it is 
and why it is useful. Naively, one would want H to be a 
true black hole event horizon. This is problematic how- 
ever as the event horizon is a global object and requires 
knowledge about the entire spacetime as opposed to a lo- 
cal picture for microstates. A local definition of horizon 
which still captures black hole thermodynamics is that 
of an isolated horizon (see also E| for overviews 
of the isolated horizons program). An isolated horizon 
is essentially just a null hypersurface with vanishing ex- 
pansion of the null normal. 1 This is still not quite good 
enough, as the approach of Carlip requires that H be an 
initial data surface and hence spacelike. If, however, one 
takes H to be a spacelike stretched horizon, which is de- 
fined as a nearly null spacelike hypersurface that traces a 
true isolated horizon, then one can use the Hamiltonian 
formulation of general relativity. 

Assuming that H is a stretched horizon, Carlip looks at 
the algebra of diffeomorphisms on H in 2-d dilaton grav- 
ity The usual generators of translations along the surface 
and local Lorentz boosts in general relativity do not pre- 
serve the horizon constraints. If H is to be a horizon, 
however, then the generators must preserve the horizon 
constraints, which means that they should be modified 
from their usual form. The necessary modification of the 
generators leads to a generator algebra equivalent to that 
of a conformal field theory with non-zero central charge 
and a calculation of the density of states using the Cardy 
formula yields the standard result for the entropy of a 



Isolated horizons also have some other properties, one of which 
is that all the field equations are satisfied on the horizon. In 
the approach presented here, the classical matter field equations 
are deliberately not satisfied on H. Therefore H is not a true 
isolated horizon, however the geometric characteristics are the 
same. 
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black hole. 

This result is intriguing as it tells us about the den- 
sity of states from the classical symmetries of the ac- 
tion, rather than some specific quantum gravity model. 
However, it also suffers some limitations. The primary 
shortcoming is that the constraint approach in the Hamil- 
tonian formulation cannot handle dynamics (as the ini- 
tial data surface is the horizon) and therefore cannot 
yield information about black hole evolution. In contrast, 
the first law of black hole thermodynamics is specifically 
about dynamics. In this paper, I show how the imposi- 
tion of horizon constraints for a dynamical Schwarzschild 
black hole generates a matter energy-momentum flux at 
the horizon that agrees with the first law of black hole 
thermodynamics. 

Intuitively how does this happen? Let us step back and 
consider more generally what happens to the structure 
of general relativity when we require H to be a horizon. 
Relativity is often labelled as "diffeomorphism invariant" 
and, of course, it is. However, special relativity is also dif- 
feomorphism invariant (as is almost any reasonable the- 
ory of physics that we have) since it produces the same 
observable results no matter what coordinate system one 
uses to calculate in. The real way in which general rela- 
tivity is unique is that there are no background geometri- 
cal objects, unlike special relativity where the spacetime 
metric is fixed to be the Minkowski metric. The fact 
that every field in the action is dynamical has important 
consequences; the most relevant of which for this paper 
is that it leads to the conservation of the matter stress 
tensor via Noether's theorem. 

The constraint approach violates this tenet of general 
relativity. The requirement that H is a horizon consti- 
tutes a fixed background geometry and hence will gener- 
ally lead to a conservation anomaly - the matter stress 
tensor cannot be conserved at the horizon. 2 If the the- 
ory is completely classical then this is a problem as the 
matter stress tensor is conserved as a consequence of the 
matter equations of motion. The only way out is for 
the matter fields to be off shell, i.e. there must be some 
purely quantum mechanical effect in the matter sector re- 
sponsible for generating matter flux at the horizon. For- 
tunately, in a black hole spacetime we have exactly such 
an effect: Hawking radiation! 

Hawking radiation from a spherically symmetric black 
hole reduces the area of the horizon according to the first 
law of black hole thermodynamics, 8ttGAQ = kAA. AQ 
is the (negative) energy flux of the current T a p£" across 
the horizon H (where £@ is the Killing field for which H 
is a Killing horizon), n is the surface gravity of the black 
hole, and A is the area of the horizon cross section. There 



2 For other e xamp les of how anomalies can be used in black hole 
physics see llCtll ll . 



is no guarantee, of course, that the flux generated from 
requiring H to be a dynamical horizon should match this 
law. However, I show that for Schwarzschild geometry 
this actually is the case up to an undetermined constant, 
and that the exact first law is recovered if this constant 
is chosen to be unity. I now turn to the derivation of this 
result. 

MODIFIED EINSTEIN EQUATIONS 

Consider a local region D of a manifold M in a co- 
ordinate system x@ = (u,v,y,z). (Eventually spherical 
symmetry will be assumed and y, z will be the angular 
directions.) I can define a hypersurface H by the condi- 
tion u = and a vector field h a such that h a \ h generates 
H . One should think of h a as the (approximate) Killing 
field for which H is to be a Killing horizon. Furthermore, 
define a covector field i a , proportional to d/dv, such that 
on H I have i ■ h = — 1. This construction is metric in- 
dependent since I have not yet imposed any restrictions 
on g a /3- I remind the reader that h a and i a are not to be 
considered dynamical fields, but rather fixed background 
geometrical objects. 

I now must decide how to impose the constraints that 
make H a horizon. A horizon emitting Hawking radia- 
tion is not quite null, but instead slightly timelike |l2| . 
Therefore the Hamiltonian formulation of Carlip is not 
directly applicable as H cannot be an initial data sur- 
face. Instead I shall work in the Lagrangian formulation 
and enforce the constraints via Lagrange multipliers in 
the action. I take H to be a timelike, instead of space- 
like, stretched horizon (a timelike membrane in the lan- 
guage of |3) • As well, I will not require the expansion 
of h a to identically vanish but instead be infinitesimally 
small, which will give an infinitesimal change in the cross 
sectional area of the horizon. Hence the constraints are 

1. $ h = h a h a +L 2 = 

2. $© = P af3 \7 a h fj -0 = 

where L, are infinitesimal on H and P a/3 is the projec- 
tion tensor P a P = g a P - h a h^/{h ■ h). 

So far I have not said anything about the surface grav- 
ity of the horizon, which also needs to be specified when 
analyzing the dynamics of a black hole. If H was a true 
null surface then L would be zero and constant along H, 
which implies that the the directional derivative of h 2 
would be of the usual form 

V a h 2 = 2nh a (1) 

where k is the surface gravity. Since h a is not quite null 
V a /i 2 is not proportional to h a . In spherical symmetry 
V Q /i 2 can be written as 

\7 a h 2 — 2nh a + pi a (2) 
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where p is a function that depends on the evolution of 
the surface. Contracting with i a gives 



i a W a h 2 



-2k- 



pi 



(3) 



If I impose the constraint i 2 = L 2 then as L — > the i 2 
term drops out of 0- So, finally, I fix the surface gravity 
k by imposing two additional constraints on H, 

1. $ 4 = « Q — L 2 = 

2. $ K = i Q V Q /i 2 + 2k = 

which gives the right value for n on the constraint sur- 
face. There is still an ambiguity in i a as the above con- 
ditions do not uniquely fix i a on D. I use this freedom 
to further choose i a such that it satisfies the geodesic 
equation, i^V pi a = 0. Note that this last choice for i a 
is not a constraint as it is not a consequence of the pres- 
ence of a horizon with surface gravity k. Therefore it 
will not be implemented with a specific constraint term 
in the action. Instead, given a solution for g a p and i a 
on H I simply choose i a to satisfy the geodesic equation 
everywhere. 

The constraints are to be imposed only on H . Using 
the Dirac delta function, defined by 



(4) 



d n xduF(u,x)S(u) = / d n xF(0,x) 



the constraints can be enforced only on H by adding 
them to the action, 

S = j \/=g~d 4 x{R + 8nG£( m) (y) (5) 

-S(u)(Xi^i + \ h $ h + A e $e + AA)] 

where C m (^) is the matter Lagrangian (as a function of 
matter fields \&) and the A's are Lagrange multipliers. 

The choice of volume element for the constraint terms 
requires some explanation. Integrating over u, the extra 
constraint terms can be rewritten as 



d 3 a 



-g(\i$i + X h $ h + A e $e + A K $, 



(6) 
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At first glance, this term seems wrong as the integration 
over H naturally involves the determinant of the induced 
metric q a p on H, rather than \J—g. However, let us con- 
sider more carefully the effect on, for example, the sur- 
face gravity constraint if the volume element was chosen 
to be \J-q_. On H , the equation of motion derived from 
varying A K is y/—q$> K = which enforces <I> K = only if 
y/—q 7^ 0. In the limit as L — ► 0, H becomes null and 
yj— q — ► as well. In this limit the constraint $ K is com- 
pletely lost as the A K equation of motion is satisfied for 
any value of $ K . Therefore the natural volume clement 
on H is not the appropriate quantity to use for enforcing 
constraints. Integrating the constraints with respect to 



the four volume element y/—g does, however, enforce the 
constraints as \J—g is non- vanishing everywhere (as long 
as the metric is non-singular). 

Varying JSJ with respect to the Lagrange multipliers 
enforces the constraints. Varying (JSJ with respect to g a ^ 
and using the identity h a V a 5(u) = (since h a is a gen- 
erator of H) yields the modified Einstein equation 



C aP = 8ttGT^\ 



C a p, the constraint "stress tensor", is given by 



C a p, = 8{u)[Xii a ip - Xhhahp + X K i^ p)h 2 

l -hahph 1 k 1 )] 



(7) 



(8) 



h 2 



1. 



+h a h V 1 {X K 8{u)F) + S(u)-V 7 (X e h^Pap) 

where h a = h^V ph a . Note that in © I have imposed 
the constraints and dropped terms that vanish (such as 
those from the variation of y/—g) when the constraints 
are satisfied. C a p can be further simplified by noting 
that in the limit of L — > 0, h a is the generator of a null 
hypersurface and so h a = fh a where / is some scalar 
function. In this limit the Ae terms in the second line of 
© cancel, leaving the simpler expression 



Cap = S(u)[Xii a ip - X h h a hp + X K i^p)h 2 
+iv 7 (A e ^ 7 P Q /3)] + h a h V y (X K 5(u)P). 



(9) 



I now discuss how the presence of C a p in the field equa- 
tions affects the conservation of the matter stress tensor. 



GENERATED MATTER FLUX 

Let us consider a simple black hole geometry, that of 
a Schwarzschild black hole. Hence from this point on 
spherical symmetry is assumed and all quantities are in- 
dependent of y,z. G a p vanishes and so contracting (JTJ 
with h 13 and taking the divergence yields 



V a (C af3 hP) = 87rGV Q (T^ l) ^). 



(10) 



I now integrate (|1(J|) over an infinitesimally wide region 
R that contains a portion of H (Fig. ^) . Since the terms 
in l|10(l are total divergences Stokes' theorem gives 



Si? 



^n a C a ph? = 8nG f V7^™V (11) 



<>FI 



where n a is the normal to the boundary dR of R and 
«/7 is determinant of the induced metric j a p on dR. I 
choose dR such that i a generates dR across H . Since 
i a is proportional to d/dv and I have the constraints 
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FIG. 1: Horizon geometry and the surface 7?. Each point 
represents a 2-sphere in (y,z). 

i 2 ,/i 2 w 0, i ■ h = — 1, i a generates translations in the 
u direction. I will therefore be integrating over u and 
so can easily evaluate the integral over 5(u). The right 
hand side of is simply 8irG times the total matter 
flux AQ emanating from the horizon through the sur- 
face dR. Without the constraints this would be zero, 
corresponding to no anomalous flux for the approximate 
Killing field h a . 

I now evaluate the left hand side of (flip. C a phP is 

C^bP = 5{u) [-X h h a h 2 -\i a (12) 

+\K(i a h?V ph 2 - V a h 2 ) - ^XePaphP] 

+h a h 2 V 1 (X K 5( u y) 

where I have applied the Liebniz rule to put the Ae term 
in the above form. In the limit as h a is null, the Ae term 
vanishes. The presence of the delta functions implies 
that there is only a contribution to the integral when 
dR crosses H (the two dimensional surfaces H- , H+ in 
Figure^. Therefore I can safely use the identity n-i — 0. 
Contracting with n a and applying n ■ i = 0, I have 

n a C a ph = -5(u)\ h h 2 (n ■ h) (13) 

--5{u)\ K n a V a h 2 + (n ■ h)h 2 \/ , 7 (5(u)A„t r ). 

I now rewrite the last term in (| 1 31) to move the delta 
function outside the derivative. Using the identity 5@ — 
la ~ n a n ^ the last term can be rewritten as 

(n ■ h)h 2 V • 1 {8{u)\ K f) (14) 
= (n ■ h)h 2 [D y (5(u)\ K F) - 5(u)X K n a n V p i a ] 

where D a = 7^V^. The last term in (|14|l vanishes in 
the limit L — * as n a w ±i a /Vi? and (as chosen previ- 
ously) i a satisfies the geodesic equation. The plus/minus 
is present as n a changes orientation from H + (where n" 



is approximately parallel to i a ) to (where the two 
vectors are approximately anti-parallel). The D 7 term 
in ITU can be integrated over dR by parts and becomes 

-<5(u)A K ^£> 7 ((n-/i)/i 2 ). (15) 

There is no boundary term due to the presence of the 
delta function. The entire integral in 1)11(1 is then pro- 
portional to S(u). Integrating over the delta function 
(|1 If) becomes 

8ttGAQ = - / ^j[\ h h 2 {n ■ h) - \\ K n a V a h 2 (16) 

-X K P 'D 7 ((n ■ h)h 2 )} 

where the notation f H H signifies that the result is the 
sum of the integrals over H + and H- . Again replacing n a 
by ±i Q / Vv 2 one sees that the A^ term vanishes as h 2 — » 0. 
Applying the constraint <i> K , h 2 — > 0, and (again) the fact 
that i a satisfies the geodesic equation the last two terms 
combine and I have 

8nGAQ= f ^\*-[ V7A K ^= (17) 
Jh + v« Jh- v« 

Now note that in this coordinate system ^/j is equal to 
Vi 2 , where 72 is the induced metric on H+ , H- . The 
factors of Vi 2 cancel and I am finally left with 

8-kGAQ = k{( - / Vt^Ak). (18) 

If X K varies between H + and H- then there is a non-zero 
AQ even if the geometry does not change. However, this 
does not correspond to the case of interest, where the flux 
is due to the time dependence of the geometry and not the 
time dependence of the Lagrange multipliers. Therefore 
A K should be chosen to be a constant. If I choose A K = 1 
then I recover that the generated matter flux satisfies the 
first law of black hole thermodynamics, 

8ttGAQ = k ([ yffi- [ ^) = kAA (19) 

JH + JH- 

which is the promised result. In summary, I have shown 
that requiring the presence of a horizon with surface grav- 
ity k by imposing constraints in the gravitational action 
also generates a (necessarily quantum mechanical) mat- 
ter flux which, at least for Schwarzschild geometry, agrees 
with the first law of black hole thermodynamics. 

DISCUSSION 

There are some open issues with the above result. The 
most pressing is obviously the value of A K . While the 
choice A K = 1 is not an exceptional or fine tuned value, 
there is no a priori reason why this choice is the correct 
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one. The difficulty in determining A K stems from the fact 
that I am using the variation of the horizon to determine 
the generated matter flux without reference to any mat- 
ter equations of motion. If the (quantum mechanical) 
matter response was known then it is possible that the 
value of X K could be fixed by applying the constraint. 
However, without the matter field equations this method 
of determining A K is not applicable. 

Another open issue is the connection with Carlip's 
method for calculating black hole entropy. One would 
hope that a connection can be made since the underly- 
ing approach, applying constraints to force a horizon, is 
the same. If a connection can be made and in the frame- 
work presented here the entropy can be shown to be A/ A, 
then the first law is a true thermodynamic equation since 
in this derivation it relates the outgoing energy flux to 
the number of microstates of the full quantum theory. 
There are a few obstacles that still need to be overcome 
in order to relate the two approaches. Since the hori- 
zon surface is a timelike hypersurface here (as it must 
be to emit Hawking radiation), it is somewhat unclear 
how to proceed with the Hamiltonian framework. Fur- 
thermore, the constraints on the surface gravity do not 
appear in the formalism of Carlip, so there is no guaran- 
tee that the algebra of constraints is identical. If these 
technical problems can be overcome, then the constraint 
program might yield not only the number of microstates 
of a black hole, but also more directly tie Hawking radi- 
ation and the laws of black hole thermodynamics to the 
microscopic theory of quantum gravity. 
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